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We develop a classical theory of electron confinement in conducting nanoparticles. The theory is 
used to compute the nonlinear optical response of the nanoparticle to a harmonic external field. 



Fundamental and applied research in the area of plas- 
monic nanodevices and of engineered materials con- 
structed from plasmonic nanoparticles is at the center 
of modern optics [l|, Most theoretical approaches 
to such systems are based on the classical electrody- 
namics of continuous media. However, in the case of 
nanoparticles whose linear dimensions are not dramati- 
cally larger than the atomic scale (a factor of ~ 10 is 
typical), finite-size and quantum corrections are of im- 
portance 0. Previously, the above two terms have been 
used interchangeably and it is accepted that the small- 
size effects are quantum mechanical in origin. Thus, 
they can be understood by considering discrete electron 
states in a nanoparticle [j, Q or reduction of interband 
screening and electron spill-over near the nanoparticle 
surface [a,0|- Li this Letter, we demonstrate that there 
is an additional, purely classical mechanism that leads 
to finite-size effects and, in particular, to nonlinearity 
of the electromagnetic response of conducting nanopar- 
ticles. We refer to this mechanisms as the classical con- 
finement effect. 

In this work, we develop a theory of classical con- 
finement of electrons in conducting subwavelength-sized 
nanoparticles and derive nonlinear polarizabilities in ID 
(slab) and 3D (sphere) geometries. Our theory is non- 
perturbative and fully accounts for electron-electron in- 
teractions within the accuracy of the quasistatic approx- 
imation. The predicted effects differ from other optical 
nonlinearities, most importantly, by the unusual depen- 
dence of the nonlinear response on the intensity of the 
incident laser beam. 

Optical nonlinearities in noble metal nanoparticles in 
different phases (such as dilute solutions, random com- 
posites and fractal aggregates) have been extensively in- 
vestigated (see Refs. @, Q and references therein). Most 
experimental studies were focused on the influence of 
the composite geometry on the nonlinear optical re- 
sponse [lfl, LLlI, 02, LUl H] • Measurements of the nonlin- 
ear response as a function of the intensity of the incident 
beam have been largely limited to nonlinear refraction 
and absorption. We note that the experimental data of 
Ref. [l2j for the dependence of the nonlinear correction to 
the absorption coefficient on the incident beam intensity 
are poorly explained and more in line with the theory 
developed below. 

On the theoretical side, the size-dependent nonlinear 



susceptibilities of a conducting nanoparticle were calcu- 
lated in Refs. [HQ from the first principles. A conduct- 
ing nanosphere was modeled as a degenerate electron gas 
placed in an infinitely- high confining potential and sub- 
jected to a time-harmonic external electric field. In this 
model, the optical nonlinearity is a consequence of the 
saturation of optical transitions between discrete elec- 
tronic states, however, the confinement effect which we 
discuss in this Letter is not taken into account. Thus, 
for example, in Ref. [5[ , a Hamiltonian was used with the 
interaction term V — er • E where E = /E cxt , E cxt being 
the external field and / = [1 + (e — 1)/3] _1 the Lorentz 
factor. In other words, it was assumed that the electrons 
move in an electric field that is obtained from the solu- 
tion to the Laplace equation which does not account for 
finite-size effects. In a more recent theoretical work @,0] , 
a finite-size effect which is conceptually related to the 
one described in this Letter but occurs in a very differ- 
ent physical setting was considered for the purpose of 
calculating the enhancement factor of Raman scattering 
by a molecule adsorbed on the surface a nanoparticle. 
However, polarization of the nanoparticle itself was not 
addressed in @, 0|- In contrast to Refs. 0, H, H, 0|, we 
present here an analytical and fully self-consistent theory 
of the nonlinear optical response of conducting nanopar- 
ticles. 

We begin by noting that the classical electrodynamic 
theory of conductors is based on the implicit assumption 
that the volume density of free charge is infinite. More 
specifically, the atomic lattice is assumed to be rigid and 
to carry a uniform positive volume charge while free elec- 
trons form a negatively-charged compressible fluid. If we 
apply an external field E cxt , the two volume charges shift 
with respect to each other by a distance 6 which results 
in the formation of a surface charge with density a. From 
the linearity of Laplace equation, it follows immediately 
that a oc 6 oc E cxt - However, the volume charge densities 
are assumed to be so large that, irrespective of the magni- 
tude of the external field, the shift 6 is much smaller than 
any other physical scale in the problem. This assumption 
is exceedingly accurate for macroscopic conductors. But 
in nanoparticles, a nonzero value of S can result in ex- 
perimentally observable nonlinearity. 

As a simple qualitative illustration, consider a conduct- 
ing sphere of radius a in an external time-varying field 
i?ext = Eq exp(— iojt) polarized along the z-axis. Here 
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and thereafter, we denote the complex representation of 
physical observables by a tilde; the corresponding real 
quantities are obtained by adding a complex conjugate, 
i.e., E = E + E* . Within quasistatics, the conventional 
result is that the sphere acquires a surface charge den- 
sity a = (3Eo/Aw)[(e - l)/(e + 2)] cos0exp(-iwi), where 
9 is the polar angle and e = e(uj) is the the permit- 
tivity of the conductor. At low frequencies, e — > ioo 
and a — (3Eq/4it) cos 9 exp(— iuit). We can use this ex- 
pression to compute the dipole moment d of the sphere, 
to obtain d = a 3 Eoexp(—iojt). Now let Eq > and 
consider such values of the time t and the angle 9 that 
cos 9 exp(— iojt) = 1. Then the real- valued surface charge 
density is a — 2a — 3Eq/2tt > 0. Since the lattice is 
viewed as rigid, the above surface charge is produced 
by a layer of depth h = 3Eq/2ttp, where p = Ze/£ 3 is 
the volume charge density of the lattice, Z is the num- 
ber of free electrons per atom, e is the electron charge, 
£ is the lattice constant, we have assumed a cubic lat- 
tice and neglected surface roughness. We thus find that 
h/a = (3/4Tr)(E /E at )(l/a), where E at = Ze/£ 2 is the 
atomic field. Once we take into account the nonzero 
value of h, the linear formula d = a 3 Eq exp(— iwi) is no 
longer valid. We emphasize that the above considerations 
applies only to the positive charge. In the classical con- 
finement model, a layer of negative surface charge with 
negligibly small depth can still be formed since the con- 
duction electrons are viewed as a compressible fluid. 

Now consider the magnitude of the effect described 
above. The theory has two small parameters: Eq/E^ 
and I /a . The first parameter is typical in nonlinear op- 
tics [15|]. The second small parameter, £/a, is negligibly 
small for macroscopic spheres. However, for a ~ 5 nm 
and £ ~ 0.5 nm (silver), the ratio is ~ 1/10. Further, we 
will see below that the nonlinearity in the classical con- 
finement model is manifested to the second order in En, 
while in the quantum- mechanical theory of Refs. 0, Q, 
the first nonlinear correction is in the third order. We 
thus conclude that the nonlinearity described in this Let- 
ter is of the same order or larger than other typical optical 
nonlinearities. 

It can be argued on very general grounds that the 
mathematical properties of the nonlinear polarizability 
in the classical confinement model are markedly different 
from that in the quantum-mechanical model where non- 
linearity is a consequence of the saturation of electronic 
transitions. To understand this, consider the polarization 
of a two-level quantum system with energy gap Hu} 2 i and 
relaxation constants T mn which is pumped with the rates 
Q rn (m, n — 1, 2; the terms T mn describe relaxation rates 
for the corresponding elements of the density matrix g; in 
a closed system Qi — T 22 g 22 ) and is excited by a time- 
harmonic field E cxt = Eo exp(— icut). The expectation 
value of the dipole moment is (d(t)) = di272.2i cxp(— iut), 
where di2 is the dipole moment matrix element and 1Z 2 i 
is the slowly-varying amplitude of the off-diagonal ele- 
ment of the density matrix. In the rotating-wave ap- 
proximation, the latter is given by 



v G Q 2 /T 22 - Qi/Tn 

21 n - iv 21 1 + 2/t(i/r 22 + i/r u ) ' 1 ' 

n = uj 21 -u;, G = -dai • E , « = . (2) 

n 2 + t 21 

If the incident wave is linearly polarized [l6|, we have 
|G| 2 = |d 2 i| 2 i?o- Thus, the expectation of the dipole mo- 
ment, as well as the polarizability tensor, are meromor- 
phic functions of the field amplitude Eq. In particular, 
(d(t)) can be expanded into a Taylor series in Eq. Due 
to the spherical symmetry of the problem, this expansion 
contains only odd powers of the field. As is the case for 
any meromorphic function whose Taylor expansion does 
not contain a zeroth-order term, the dipole moment van- 
ishes in the limit Eq — > oo. Physically, this corresponds 
to saturation of the transition. Indeed, the second factor 
in the right-hand side of Eq. ([1]) describes field-dependent 
population inversion which, in the limit of a strong in- 
cident field, approaches zero. The saturation process is 
much more complicated in the case of a degenerate elec- 
tron gas in a spherical potential well which was consid- 
ered in Refs. [j, H|, yet the analytical properties of the 
solutions are essentially the same. 

In the case of classical confinement, the situation is 
substantially different. Here the dipole moment does not 
vanish in the limit of strong external fields but rather 
oscillates with an amplitude that approaches a constant 
limit. Consider again a conducting sphere in a time- 
harmonic external field. We will ignore the possibility 
of ionization, i.e., similarly to Refs. 0, @, we will as- 
sume the existence of an infinite potential barrier that 
prevents electrons from leaving the system. Then, for a 
very strong external field, the amplitude of the oscillat- 
ing dipole moment is d max = Nea, N being the number 
of free electrons in the system. Here d ma x is the dipole 
moment which is obtained when all conducting electrons 
are concentrated at a point near the sphere's surface. 
Although the limit we have just considered is arguably 
unphysical, its mathematical implications are important. 
Indeed, there is no such analytic function of Eq that is 
zero at Eq = and approaches a finite limit at infinity. 
We conclude that the function d(Eo) is not analytic and 
can contain, for example, a term oc Sol-^ol which, al- 
though it is second order in Eq, docs not need to vanish 
in centrosymmetric systems. We will see below that this 
is, indeed, the case. 

We now proceed with more detailed calculations. 
As a first step, consider the one-dimensional problem 
schematically illustrated in Fig. la. Here an external 
field -Ecxt = —Eq is directed perpendicularly to a slab 
of thickness L. If Eq > 0, the free charge is distributed 
inside the slab as follows. The surface which is opposite 
the field direction (as shown in the figure) acquires a 
negative surface charge a 2 — —Eq/Att. Near the other 
surface, a positively-charged layer of depth h = Eq/Attp 
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is formed, where p = Ze/l 3 is the background positive 
volume charge density. We thus see that the slab is sep- 
arated into two regions. The first region is characterized 
by zero conductivity due to the absence of free carriers 
(this region is dashed in Fig. la) while the second region 
is conducting. In static equilibrium, the local field in 
the second region must be zero, while there is no such 
requirement for the first region. We note that the field 
inside the conducting region which is produced by the 
positively charged layer is 2-ko\ where <j\ — E /4ir, 
the same as would be produced by a surface charge 
of density a%. Thus, the depolarizing field inside the 
conducting layer is -Bdep = 2it(ai — era) = Eq and the 
total local field E\ oc = E cxt + Ed cp is zero. 
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FIG. 1. Illustrating the geometry of the ID (a) and 3D 
(b) problems. 

Next, consider a time- varying external field E cxt (t) = 
Eq exp(-iujt). The positively-charged layer can now ap- 
pear on either side of the slab. It is convenient to in- 
troduce the following notation: let the total charge per 
unit area which accumulates near the left face of the slab 
be u\ and the similar quantity near the right face be 



(72- Here we do not distinguish between a true surface 
charge and a surface layer of finite depth. It follows 
from charge conservation that <j\ = —02- We then can 
write o\ = —a, 02 = +c, where a can be either positive 
or negative, depending on the phase of the oscillations. 
The depolarizing field in the conducting region is given 
at any time by E^ep = — 47rcr and the total local field by 
-Bloc = E ext + E^cp- Then the equations of motion for a 
negative test charge inside the conducting region can be 
written as 



m[v 



-e[E ext (t) - 47rcr] 



-pv 



(3) 



Here v,m and e are the electron velocity, mass and 
charge, respectively, and 7 is a phenomenological friction 
term; the first equation in © is Newton's second law 
while the second is the continuity equation. The oscilla- 
tory solution to the above system of ordinary differential 
equations is 



4Tra(t) = ujlE e- iuJt 



H 



■ to 



(4) 



where lo p = ^jA^ep/m is the plasma frequency. Wc 
can further compute the local field i?i oc and the current 
j = —pv inside the conducting region and verify that the 
ratio j{t)/E\ oc {t) yields the classical Drude conductivity 
iuip/4Tr(ui + ij). 

So far, the results appear to be conventional. Finite- 
size effects and the nonlinearity of the optical response 
become apparent when we compute the dipole moment 
per unit area of the slab, V — Ad/ AS. A straightfor- 
ward calculation yields V(t) = cr(t)[L — h/2]. Thus, 
the effective width of the slab is reduced by h/2 where 
h = \<j\/ p (note that h is related to the absolute 
value of the real-valued quantity a). We now find that 
V{t) = <j{t)[L - (l/2p)|<r(t)|]. Wc further note that 
the field amplitude E can always be chosen to be real 
so that a(t) — EqA(ui) cos(cut — if), where A(lo) = 

uj 2 p /2it^J {uj^ - lu 2 ) 2 + (7W) 2 and tarn^ 

We then write the final result as 



Juj/{uj 2 



V(t) = aE cos(ujt — tp) 1 + (3\E cos(wi — (p) 



(5) 



where a = LA(u), f3 = -A(u)/2pL = -A{uj){l/2L)E~ t 1 
and we have introduced again the microscopic quantities 
£ = {Ze/p) 1 / 3 and E at = Ze/l 2 . 

Note that ([5]) is not an expansion but is exact as long 
as h < L or, equivalently, max(|er|) < pL. The impor- 
tant feature of the obtained solution is that "P(-Eo) is 
not an analytic function and can not be expanded into 
a Taylor series. This mathematical property is closely 
related to the existence of a finite limit lim^^co V(Eq) 
(in the saturation model of Refs. 0, [f| this limit is zero). 
The non-analyticity of the optical response has important 
consequences for the generation of frequency harmonics 
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as is shown below for the more physically important 3D 
case. 

We now consider the problem of a three-dimensional 
conducting sphere. It turns out that accounting for clas- 
sical confinement in this case leads to a formidable math- 
ematical problem. We will use, however, certain physical 
insights that will allow us to obtain a nonperturbative 
analytical theory. Consider a conducting sphere of ra- 
dius a, a constant background positive volume charge 
p and a free charge whose integral over the sphere vol- 
ume is — 47ra 3 p/3; an infinitely high spherical potential 
barrier that prevents ionization is assumed. We seek to 
find the dipole moment of the sphere in an external field 
E oxt = zEq enp(-iujt) , where z is a unit vector pointing 
in the direction of the z-axis. As before, we recognize 
that the sphere is separated at all times into two regions: 
one region has no free carriers and is nonconducting while 
the other region has a constant nonzero conductivity; this 
conducting region is doubly dashed in Fig. lb. 

The first physical observation that we make is that in 
the quasistatic problem with a time-harmonic external 
field, the motion of charges is such that, at any time t, 
both the volume and the surface charge densities corre- 
spond to a static equilibrium which is obtained for an ex- 
ternal field A(w)E ex t(i') which is taken at a different time 
t' and multiplied by a frequency-dependent real-valued 
factor A(ui). Thus, the system goes through states of 
static equilibrium which are phase-shifted with respect to 
the external field. Mathematically, this statement follows 
from the linearity of the equations of motion. In static 
equilibrium, the electric field in the conducting region is 
zero. In the dynamic problem, the latter is nonzero but 
proportional to the difference E ext (t) — E ex t(*')- We as- 
sume here that E oxt is spatially homogeneous over the 
volume of the sphere. The only motion of the free charge 
inside the conducting region that is consistent with this 
condition is one-dimensional motion along the z-axis. 
From this, we find that the surface that separates the 
conducting and non-conducting regions is a sphere. The 
center of this sphere is denoted by O' and is shifted from 
the center of the original sphere by a distance S along 
the z-axis, where S can be both positive and negative 
(see illustration in Fig. lb). We thus can characterize 
the dynamics of the system by a single scalar parameter 
5. 

The second observation will allow us to find the depo- 
larizing field inside the conducting region. As we have 
argued above, the total electric field inside that region 
is spatially homogeneous and directed along the z-axis. 
The external field does satisfy this condition and so must 
the depolarizing field Edcp- The latter is a superposition 
of the field produced by a positively charged meniscus 
and the negative surface charge a which we have not yet 
determined. We notice, however, that a field with the 
required properties is created by two oppositely charged 
menisci of the same shape which are shown in Fig. lb 
as single-dashed regions. Indeed, the field produced by 
the two menisci is the same as the field of two oppo- 



sitely charged spheres shifted with respect to each other 
by a distance S. Inside the conducting region, this field is 
given by Ed cp = AnpSz/i. With the understanding that 
this field is created by the positively charged meniscus 
and by a yet unknown negative surface charge a whose 
field in the conducting region is the same as that of the 
hypothetical negatively-charged meniscus, we write the 
equation of motion as 

m (d + Y<5) = -e [z • E ext (t) + iirpS/3] . (6) 
The oscillatory solution to ^ is 



5{t) = -(e/m)E e~ 



\ujp — uj 2 — ijuj] 



(7) 



where ujp = w p /v3 is the Frohlich frequency. In the con- 
ventional approach, the dipole moment of the sphere is 
computed as d — — Airpo?5/3. Evaluation of this expres- 
sion leads to the linear polarizability a — a?(e — l)/(e+2) 
with e = 1 — ujp/uj(uj + ij). We, however, intend to 
take into account the presence of the meniscus and the 
fact that the surface charge density can deviate from 
the usual oc cos 9 dependence. To this end, we write 
d z = z • d = p f v zd 3 r + f g zad 2 r. Here the first inte- 
gral is over the volume of the positively charged meniscus 
while the second term is the contribution of the negative 
surface charge. After tedious but straightforward inte- 
gration, we obtain 



2tt 



pS 
3 



11 1 4 16 



a 3 ai 



(8) 



where o\ — J cr(cos ff) cos 6d cos 9. We thus see that 
knowledge of the first moment of the surface charge den- 
sity suffices for the purpose of computing the dipole mo- 
ment. While finding the whole function ^(cos^) is a dif- 
ficult task, <7i can be found by considering the depolar- 
izing potential in the vicinity of the sphere origin, O. 
On one hand, we know that the potential is given by 
the formula (j) — —E^cpZ with Ed ep — Airp5/3. On the 
other hand, we can write the same potential as an integral 
over the meniscus and the surface of the sphere and ex- 
pand the resultant formula into scalar spherical harmon- 
ics r Yi m (i • z). A straightforward calculation of the term 
I = 1, m = yields <f> 10 = 2tt[cti - (p<5/3)(l + (3|5|/8a))]z. 
[The condition that all higher terms in the above ex- 
pansion vanish provides an infinite set of equations for 
the higher moments of cr (cos 6).] From the equality 
<?Wp = 0io j h follows that 

en = -(p8/3) [1 - (3\5\/8a)] . (9) 
We then substitute the above result into (El) to find 



d z = 



AirpSa 3 



, 9 W 1 a 2 ]*] 
16 a 32 a 3 



(10) 
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where 8 = 8+6* and 8 is given by ||7J). We emphasize that 
the above formula is not an expansion: it is exact as long 
as \S\ < a, the latter condition having been implicitly 
used for computing a±. 

We thus find that the polarizability of the sphere has 
a mathematical structure which is similar to that of 
the one-dimensional slab, except for the additional term 
oc 8 3 \8\ which appears in (|10p . For most practical pur- 
poses, this term can be neglected since it contains an ad- 
ditional small factor (E$l/ E at a) 2 compared to the term 
tx 8\5\. Both terms, however, describe the generation of 
odd-order frequency harmonics, i.e., oscillations at the 
frequencies u> n — (2n + l)u>, n — 0,1,2,.... Indeed, 
let Eo be purely real, so that S\8\ oc cos(r)| cos(r)| and 
8 3 \8\ oc cos 3 (r)| cos(r)|, where r = ut — ip, and consider 
the following Fourier series: 

r \\ i \\ ^ 8sin(7rfc/2) 
cos(r)| cos(r)| - 7rfc(4 _ fc 2) cos ( fcr ) ( n ) 

,, , . , 96sin(7rfc/2) . 

cos3(,)|cos(r)Hg , fc ( 64 ,y + ^) C ° S(fcT) (12) 

In both cases the Fourier coefficients are nonzero only if 
k = 2n+l (the k = 2 Fourier coefficient in Eq. (fTTj) is zero 
despite the existing uncertainty). Note that even though 
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the first nonlinear correction is second order in Eq, there 
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In summary we have described novel nonlinear opti- 
cal effects in conducting nanoparticles and developed a 
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ducting nanosphere (Eqs. (|T|. (fTU)) -(fT2" )) h Experimentally 
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